from Excitation_Carbon import *

cs = CS_Excitation(species_name = "C")

id = 0

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 1.26,
    A   = 1.277,
    B   = 1.304e1,
    C   = -3.314e1,
    D   = 6.357e1,
    E   = -4.621e1,
    F   = 4.0e-2,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p21D")

#the result data is not same as the ref
id = id + 1
cs.add_cross_section(
    id = id,
    Vif = 2.68,
    A   = 1.598e-1,
    B   = 2.836,
    C   = -9.813,
    D   = -1.012e1,
    E   = 1.233e-1,
    F   = 0.0,
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p21S")


id = id + 1
cs.add_cross_section(
    id = id,
    Vif = 4.18,
    A   = -2.302,
    B   = 8.923,
    C   = -2.873e1,
    D   = 4.555e1,
    E   = -2.431e1,
    F   = 1.323e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s2p35S")

        


id = id + 1
cs.add_cross_section(
    id = id,
    Vif = 7.48,
    A   = -6.117,
    B   = 2.989,
    C   = 4.148,
    D   = 0.0,
    E   = 1.222e1,
    F   = 0.0,
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3s3P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 7.68,
    A   = 2.423,
    B   = 7.590e-1,
    C   = -6.486,
    D   = 3.143,
    E   = -2.074e1,
    F   = 6.744e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p3s1P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 7.94,
    A   = 1.307e1,
    B   = -3.052e1,
    C   = 1.835e1,
    D   = 0.0,
    E   = 1.001e1,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s2p33D")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 8.53,
    A   = 1.909,
    B   = -6.101e-1,
    C   = -2.336,
    D   = -3.093e2,
    E   = 7.450e2,
    F   = 1.491,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p3p1P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 8.64,
    A   = 9.753e-1,
    B   = 1.908,
    C   = 2.161,
    D   = -4.891,
    E   = 0.0,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3p3D")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 8.77,
    A   = 1.886e-3,
    B   = 2.263e-1,
    C   = 8.836e-1,
    D   = -1.031,
    E   = 0.0,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3p3S")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 8.85,
    A   = 7.728,
    B   = 8.661e-1,
    C   = -1.640e1,
    D   = 7.754,
    E   = 0.0,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3p3P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.00,
    A   = -1.389,
    B   = 1.062,
    C   = -3.742e-1,
    D   = 4.786,
    E   = -4.126,
    F   = 2.606e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p3p1D")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.20,
    A   = 3.567e-2,
    B   = 5.081e-1,
    C   = -7.062e-1,
    D   = 1.463,
    E   = -2.569,
    F   = 4.101e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p3p1S")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.33,
    A   = 1.146e1,
    B   = -2.477e1,
    C   = 1.281e1,
    D   = 0.0,
    E   = 6.328,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s2p33P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.63,
    A   = 5.747e-2,
    B   = -1.444e-1,
    C   = 1.793,
    D   = -1.719,
    E   = 0.0,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3d1D")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.69,
    A   = -3.067,
    B   = 2.352,
    C   = -3.127,
    D   = -9.123,
    E   = 1.656e1,
    F   = 1.476e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p3d3F")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.71,
    A   = 4.081,
    B   = -9.274,
    C   = 5.423,
    D   = 0.0,
    E   = 4.970,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3d3D")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.73,
    A   = 1.183,
    B   = -8.798e-1,
    C   = 7.839,
    D   = -1.718e1,
    E   = -8.149e-1,
    F   = 5.66e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s22p3d1F")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 9.83,
    A   = 1.247,
    B   = -5.247,
    C   = 3.983,
    D   = 0.0,
    E   = 3.876,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s22p3d3P")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 12.2,
    A   = -3.820,
    B   = 2.082,
    C   = 2.801,
    D   = -9.506,
    E   = 1.552e1,
    F   = 2.937e-1,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s2p31D")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 13.1,
    A   = 8.801e-1,
    B   = -9.630,
    C   = 8.179,
    D   = 0.0,
    E   = 1.199e1,
    F   = 0.0, 
    omega = 9.0,
    equation_number = 6,
    reaction_name = "Excitation_C_2s22p23P-2s2p33S")

id = id + 1
cs.add_cross_section(
    id = id, 
    Vif = 14.9,
    A   = 2.975,
    B   = -4.065,
    C   = -2.943e-1,
    D   = 7.361e1,
    E   = -3.359e2,
    F   = 1.40,
    omega = 9.0,
    equation_number = 7,
    reaction_name = "Excitation_C_2s22p23P-2s2p31P")

cs.plot_cs()
cs.save_cs()
cs.cal_total_cs()

